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Abstract 

The block canonical coset decomposition is developed as a universal algorithmic tool to synthesize 
n-qubit quantum gates out of experimentally realizable NMR elements. The two-, three-, and four- 
qubit quantum Fourier transformations are worked out as examples. The proposed decomposition 
bridges the state of the art numerical analysis with NMR quantum gate synthesis. 
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Introduction. According to the quantum computing circuit model, a quantum gate is 
a unitary operator designed to perform a predefined operation. The primary aim is to 
synthesize unitary operators suitable for experimental implementation. Currently, quantum 
gate synthesis is achieved either by factorization methods, as in the strategy "divide and 
conquer," or by optimal control theory [IH3]. 

The factorization methods' development was initiated by utilizing the KAK Cartan de- 
composition to attain minimal time two-qubit operations [3H6] with subsequent extensions 
to multi-qubit gates [51 U\- However, the Cartan decomposition, coming from Lie group 
theory, was not constructive and actual applications required the independent development 
of numerical routines [8j |9]. Inspired by numerical analysis, alternative approaches, built 
upon the cosine-sine matrix factorization , were employed to constructively decom- 

pose an arbitrary quantum operation in terms of single-qubit and CNOT gates as well as to 
estimate the complexity of such implementations. Since CNOT gates are not most conve- 
nient in NMR setups [T5] , these developments did not find experimental applications. While 
NMR two-qubit synthesis has been thoroughly explored (U HH ELT] , multi-qubit factorization 
remains largely an open problem [T8H20] . 

In this Letter, we introduce the canonical coset decomposition as a constructive fac- 
torization technique to build quantum gates out of standard single- and two-qubit NMR 
operations. To illustrate this method's ability, we present the explicit factorizations of the 
two-, three-, and four-qubit quantum Fourier transform, which is the basic component of 
many quantum computing algorithms. 

Canonical Coset Decomposition is a general group-theoretic tool originally devised to 
study the geometric structure of unitary matrices [2TJ |22] . In particular, this method allowed 
to obtain analytic expressions for the measure of unitary operators [211 [22] as well as the 
measure [23J and metric [23] of mixed quantum states. 

The canonical coset decomposition of a unitary matrix U is 
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where Uu and [1 — f/ji^i] 1 ^ 2 are invertible. Throughout the Letter we exclusively employ 
the block canonical coset decomposition, which is a special case of Eq. Q with U ij} X, and 
Vi,2 being square matrices of the same dimension. 

The block canonical coset decomposition is a a very powerful tool because each term in 
the r.h.s. of Eq. is simpler than the original matrix U. In particular, the first term can 
be expressed for X = 
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additionally, the second term forms a subgroup. In other words, the Lie algebra C(U) of 
unitary matrices is decomposed into a block diagonal subalgebra t and its complement p 
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For n-qubit gates U(2 n ), the factorization ([TJ) reads 
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U^- 1 ) ® U(2 

where the right factor corresponds to the block diagonal subgroup. 

The canonical coset decomposition is closely related to other techniques. In particular, 
if X is diagonal then the first term in the r.h.s. of Eq. ([!]) is equal to the middle term 
of the cosine-sine decomposition [ID] . When X is a single column, the decomposition (jTJ) 
is equivalent to the Householder transformation [25], which is not only one of the most 
important operations in numerical analysis but also very useful for the generation of 
n-level unitary operators 



3 



(A) 




5 00 


0"! $ 


5 <Tl 




5 cr 2 






02 f < 




02 $ 


5 a\ 




$ (72 


T 2 (8) <T3 


(B) 




5 cr 




S> 0"! 


o"o 5 


$> (T 2 


to (8) 0-3 




03 G 


5 ctq 


cr 3 (g> cri 


cr 3 (g> cr 2 


0-3 ® 0-3 



TABLE I: Basis elements for (A) the complement p and (B) the subalgebra t of the u(4) Lie 
algebra. 



Two-qubit gates. According to Eq. ^8j), the £7(4) group associated with two-qubit gates 
is decomposed as 

£7(4) 
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The reverse form of Eq. (|9| , 



£7(4) = £7(2) ® £7(2): 
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£7(2) ® £7(2) ■ 

can be achieved by using the identity £7 = (£7^)1 and performing the decomposition ([9J on 
£7t. 

The matrix log of £7(4) is spanned by the basis elements listed in Tables [JA) and[l]^B), 
where the basis of the subalgebra t is given in Table^B). Since t includes both single and two- 
qubit interactions, isolation of single-qubit operations necessities the following factorization: 

£7(2) ® £7(2) . 
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where S2 = [GiG^] 1 ^ 2 and S% = S\G\. Equation (13) accomplishes the desired isolation. 



Even thought Eq. ( 14 ) is not a single-qubit operator, it can be readily transformed into one 



by employing the standard NMR methods. Thus, making the factorization (11) suitable for 
NMR experiments. 



The first factor in the r.h.s. of Eq. ^ can be further simplified by a recursive factorization 
Q and g0|) 
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where the middle term becomes simple enough for NMR synthesis. Having initiated the 
unitary operators 
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the recursive decomposition is realized by the following algorithm: 

1. Perform the factorization (B of U and assign [/« := u{ ^ i2) , t/ {2) := 17(2) <g) C/(2); 

2. £7(2) := t/( 2 )?7(2), f/ ,(1) := V^U^V: 



3. Perform the decomposition (hoj of £/ ,(1) and assign := V[C/(2) g> U{2)]V\ W { - 
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4. WW := H^ilfW, U := W^; 

5. Return to step 1 if the convergence conditions ~ 1 and W^ 1 ' ~ 1 are not met; 



6. Finally, the right subgroup in Eq. (15) is stored in and the left subgroup - WW. 



The middle term U(i)/[U(2) <g> ?7(2)] 2 is in W^ 2 ). 



The matrix log of the left, middle, and right terms in Eq. (15) are spanned by the elements 



from Tables [TT]( A) , (TTJ^B), and [TT](C) , respectively. Since the middle term is made of the 
elements {a 2 <8> eri}i=o,i,2,3, it is wen suited for NMR synthesis. 

Importantly, the developed algorithm is applicable to an arbitrary n-qubit system after 
adjusting the definition of V according to 
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(n — 1) times 

T/ie two-qubit quantum Fourier transform F is used to illustrate the proposed algorithm. 
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TABLE II: Basis elements for the Lie algebra u(4) corresponding to (A) left, (B) middle, (C) right 



terms in Eq. ( 15 ). 



As shown in Sec. I of Ref. |30j, the recursive decomposition (15) leads to 
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with 

zlog(Fi) =0.55536 (T ® (01 + 02), 
z log(F 2 ) = - 0.392699 ct 2 ® 00, 
i log(F 3 i) = 0.55536 a ® (<7i - cx 2 ), 
i log(C/ 32 ) = O ® (-1.1781ao + 0.785398<x 3 ). 
The pictorial representation of this factorization is shown in Fig. [TJ 
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FIG. 1: The schematic representation of the two-qubit Fourier matrix factorization (18). Black 



boxes represent single-qubit 7r/4 rotations; gray boxes - taylored single-qubit operators specified 



byEq. (19) 



Three-qubit gates. According to Eq. (|8j), a three-qubit gate is decomposed as 
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where Z7(4) <g> Z7 (4) is block diagonal. The two-qubit decomposition (15) is extended to 
three-qubit gates 
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where the central factor's matrix log is spanned by {&2ij}i,j=o,i,2,3, where a^k = (J% ®Oj ®o k 



According to Eq. (17), the factorization algorithm is applicable in the three-qubit case once 



the definition of V in Eq. (16) is replaced by V ■= exp(i7T(72oo/4). The decomposition (11) 
is naturally generalized to isolate the subgroup U (4) associated with two-qubits 
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where the matrix logarithms of U(A) and U(A) ® U(A)/U(A) are linear combinations of 

Woij}i,j=o,i,2,3 and {o"3jjkj=o,i,2,3, respectively. 

The factorization plj) of the three-qubit quantum Fourier transform results in F = 
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U-JJ 2 \] 3 . Each of these factors are further decomposed in Sec. II of [30] 
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where Fj..., Wj... and Ujk- are single-qubit operators. This factorization is also depicted in 
Fig- § 

Finally, the four-qubit Fourier transform is elaborated in Sec. Ill of Ref. |30j . 

Conclusions. The block canonical coset decomposition is introduced as an algorithmic 
procedure to synthesize an arbitrary n-qubit quantum gate out of one- and two-qubit NMR 
operations. In particular, we worked out the implementations of the two-, three-, and 
four-qubit quantum Fourier transform given by a dense matrix (i.e., without zero entries), 
therefore presenting a convincing illustration of our method's capability. Moreover, a highly 
optimized numerical implementation of the developed algorithm can be carried out through 
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FIG. 2: The schematic representation of the three-qubit Fourier matrix factorization from Eqs. 
(23), (24), and (25). Black boxes represent single-qubit 7r/4 rotations; gray boxes - taylored 
single-qubit operators. 

the block Householder decomposition [31TI34"] . Thus, the current work provides an important 
bridge between the state of the art numerical analysis and NMR quantum gate synthesis. 
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